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We investigate non-standard neutrino interactions (NSIs) in the Zee-Babu model. The size of 
NSIs predicted by this model is obtained from a full scan over the parameter space, taking into 
account constraints from low-energy experiments such as searches for lepton flavor violation (LFV) 
and the requirement to obtain a viable neutrino mass matrix. The dependence on the scale of new 
physics as well as on the type of the neutrino mass hierarchy is discussed. We find that NSIs at the 
source of a future neutrino factory may be at an observable level in the v e — > v T and/or — > v T 
channels. In particular, if the doubly charged scalar of the model has a mass in reach of the LHC and 
if the neutrino mass hierarchy is inverted, a highly predictive scenario is obtained with observable 
signals at the LHC, in upcoming neutrino oscillation experiments, in LFV processes, and for NSIs 
at a neutrino factory. 



I. INTRODUCTION 

Experimental studies of neutrino oscillations have pro- 
vided us with compelling evidence that neutrinos have 
masses and lepton flavors mix. Among many possible 
mechanisms to describe the origin of neutrino masses, ra- 
diative mass generation provides an attractive method to 
obtain small neutrino masses. In such a framework, neu- 
trino masses are exactly vanishing at tree level, and are 
induced as finite radiative corrections. Typically, neu- 
trino masses are suppressed by a loop factor and pro- 
portional to firnj/Al 2 , where mg are the charged lep- 
ton masses, M is the mass scale of the new particles 
in the loop, and p is the scale of lepton-number viola- 
tion. Toghether with a modest suppression from Yukawa 
couplings, this allows sub-eV neutrino masses, while hav- 
ing new physics not too far from the electro-weak scale, 
M ~ 1 TeV, opening the possibility of collider tests of 
the neutrino mass generation mechanism. 

An economical way of radiative neutrino mass gen- 
eration is to enlarge the scalar sector of the Standard 
Model [l|, 0. In the Zee model, neutrino masses are 
obtained at one-loop level by adding a singly charged 
scalar and a second 5[/(2)l doublet to the Standard 
Model [|[ . While the simplest version of this model can- 
not accommodate current experimental data, since the 
predicted leptonic mixing angle #12 is too large (close to 
7r/4), a minor extension of the model remains a viable op- 
tion Q . Alternatively, in the Zee-Babu model [f| @, 0] , 
two 5J7(2)l singlet scalars are introduced, one singly 
and one doubly charged, and neutrino masses are gen- 
erated at two-loop level. Phcnomcnological studies of 
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this model have been performed, e.g., in Refs. (11. [oL fl0| ■ 
Through the exchange of heavy scalars, lepton flavor vi- 
olating (LFV) processes such as pi — ► 3e and pi — ► &y 
can be dramatically enhanced compared to the Standard 
Model. Furthermore, the new scalars could be accessi- 
ble for the Large Hadron Collider (LHC). In particular, 
the doubly charged Higgs may induce very clean like-sign 
bi-lcpton events. 

Besides colliders, next generation neutrino oscillation 
experiments will also help us to unveil the underlying 
physics behind neutrino masses. A neutrino factory 
will be the ultimate facility to perform precision mea- 
surments of standard neutrino oscillations as well as to 
search for non-standard neutrino properties. We take 
this as a motivation to investigate non-standard neu- 
trino interactions (NSIs) in realistic neutrino mass mod- 
els. In a given model, NSIs are typically linked to LFV 
for ch arg ed leptons, yielding too tight bounds, see, e.g., 
Refs. [U, [13, EH- For example, in the case of triplet 
scalar models (i.e. type-II seesaw), NSIs are always en- 
tangled with the interactions among four charged lep- 
tons, which suffer stringent constraints from LFV pro- 
cesses like fj, — ► 3e. There are no sizable NSIs unless 
severe fine-tuning of Major ana phases is invoked (l4l |. In 
the case of the Zee-Babu model, the situation is more 
involved, since the masses of singly and doubly charged 
Higgs in principle can be well separated and a differ- 
ent set of Yukawa couplings controls charged lepton and 
neutrino interactions with the scalars. Non-trivial pre- 
dictions for the NSI parameters emerge from the differ- 
ent combinations of Yukawa couplings responsible for the 
neutrino masses and mixing, and LFV processes. In this 
respect, we will investigate NSIs in the Zee-Babu model 
in detail. 

The remaining parts of the work are organized as fol- 
lows: In SecHH we sketch the Zee-Babu model and show 
how NSIs arc induced. In Sec. IIII1 experimental con- 
straints from low-energy observables on the model pa- 
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ramcters arc summarized, while in Sec. IIV1 we present 
the results from our numerical study of NSI parameters 
within this model. Discussion and summary follow in 
Sec. El 



II. THE ZEE-BABU MODEL AND 
NON-STANDARD NEUTRINO INTERACTIONS 

In the minimal Zce-Babu model, two 5{/(2)l singlet 
scalar fields h + and k ++ are introduced with hyper- 
charges 1 and 2, respectively. The corresponding La- 
grangian is then given by 



C = £sM + fapLj ia Ci(J2LLl3h + 

- iih~ h~ k ++ + h.c. + Vh , 
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where Ll denote left-handed lepton doublets, e are the 
right-handed charged leptons, and the scalar potential 
Vh contains the couplings among scalar fields. The 
Yukawa couplings / and g are antisymmetric and sym- 
metric, respectively. The trilincar /i term breaks lepton 
number (L) in an explicit way, and hence, one can natu- 
rally expect the dimensionful parameter fi to be reason- 
ably small, since the symmetry is enhanced in the limit 

Light neutrino masses are generated via a two-loop di- 
agram, which yields 



{m v ) ab = 16/J,facm c g* d l c dmdfbd 



(2) 



where m c are charged lepton masses and I c d is a two- 
loop integral Since the e + e" collider LEP at CERN 
indicates that the masses of charged scalars are typically 
larger than 0(100 GcV), we can neglect the masses of 
charged leptons compared to them. In this case, one 
finds 
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where M = max(mt, m h) and I(r) is a dimensionless 
function of order unity [lOj. For our numerical calcula- 
tions we use the expression given in Eq. (7) of Ref. Q. 
Using Eq. ([3|), the light neutrino mass matrix becomes 
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where the matrix D e = diag(m e , m^, m T ) contains the 
charged-lepton masses. Light neutrino masses are sup- 
pressed by the heavy scalar masses and proportional to 
the lepton-number violating parameter /x. Due to the 
antisymmetric property of /, we have det m v = 0, and 
therefore, one of the neutrinos is massless if higher-order 
corrections arc not considered. 

The heavy scalars will induce non-standard lepton in- 
teractions via tree-level diagrams as shown in Fig. [1] 
After integrating out the heavy scalars, the following 
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FIG. 1: Tree- level diagrams with the exchange of heavy 
scalars. The corresponding diagrams are responsible for (a) 
non-standard interactions of four charged lepton, and (b) non- 
standard neutrino interactions. 



dimension-6 operators are generated at tree level [1^ 
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where a Fierz transformation has been applied for the 
second step, and 



p pa - 



Upf. 



V2G 



F m\ 



VTeV/ 



(6) 



arc the canonical NSI parameters. For neutrinos prop- 
agating in normal matter, only the following NSIs are 
induced 
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Taking into account the antisymmetric property of/, one 
can find that the only relevant NSI parameters in matter 
are e™ , e™, and e™ . Furthermore, NSIs may show up at 
neutrino production in a neutrino factory, related to the 
processes fi — > eT7pv a due to e e J^. To be consistent with 
the notation in the literature, e.g., Ref. [l6[, we define 



r' r = 



ffJ,ef eT 

V2G 

JfiT fen 

\/2G 



F m 2 h 



F m\ 



(8) 



which correspond to the source effects in the — > v T 
and v e — > v r channels, respectively. By definition, the 
relation 

= -t 5 '* (q\ 

holds, since both the NSI parameters e™ T and are re- 
lated to the Yukawa couplings f efl and f eT . Let us men- 
tion that Eq. ([9]) holds in a rather general class of models, 
where NSIs are generated by dimension-6 operators [l?] ] . 

The light neutrino mass matrix can be diagonalized by 
means of a unitary transformation as 



UDU T 



(10) 
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where D = diag(mi, 7712, TO3) and U can be parametrized by using three mixing angles and two CP-violating phases 
I 
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where Cy = cos%, = sinfty (for ij = 12, 13, 23), 
and 5 is the Dirac CP-violating phase. Here only one 
Majorana phase a is involved, since one light neutrino is 
exactly massless. As a result of det/ = 0, there is an 
eigenvector vq = (f^r, — ferilen) which corresponds to 
the zero eigenvalue fva = Note that vo is also an 
eigenvector of m„, and therefore, we have 



e'" are real, whereas the phase of e s e 
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for both normal mass hierarchy (mi <C to 2 <C TI3, NH) 
and inverted mass hierarchy (ni2 > mi ^> 7713, IH). 

Equation l[T2"|) provides us with three equations, among 
which one is trivially satisfied, since one element of D is 
zero. The other two equations lead to relations between 
/ and the lepton mixing parameters. In the NH case, we 
have 
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According to the current global fit of neutrino oscillation 
experiments, the second terms in the right-hand sides of 
Eqs. I|13p and (|14p can be neglected, since they are sup- 
pressed by the small mixing angle 6*13. Then, we obtain 
the approximate relation f e)i ~ f er ~ Taking 
into account the experimental constraints from Eqs. Q18|) 
and (fT9")l below, we can roughly estimate that |/ eM | ~ 
|/e T | < 0.05 (m h /TeV) and |/, lT | < 0.1 (m ft /TeV). Com- 
pared with Eqs. (O and JSJ), the possibly important NSI 
parameter is e s eT , which mainly affects the v e — ► v r chan- 
nel. The CP-violating phases of f a p are suppressed by 
6*i3, and thus, NSIs cannot induce very distinctive CP- 
violating effects in the NH case. 

For the IH case, the two non-trivial equations are 
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Indeed, it is obvious that |/ e/i | ~ |/ eT | and |/ MT | ~ 
|/ct|si3/s23 hold. Thus, the potentially sizable NSI pa- 
rameters are e™ , e™ , e™ for neutrino propagation in 
matter, and e* T for source effects in the — > v T chan- 
nel at a neutrino factory. Equations (|15p and (jTSJ) imply 



that e* and 

given by 5. This may lead to an interesting correlation 
of CP-violation in standard oscillations and £g T -induccd 
CP-violating effects [HG1]. 



III. EXPERIMENTAL CONSTRAINTS 

At low-energy scales, stringent constraints from LFV 
processes mediated by the heavy scalars have to be in- 
cluded when we confront the model with experimental 
data. In the following, we compile the bounds given in 
Ref. [nj. 

• £~ — > £££~£j. As shown in Fig. [TJ these rare lep- 
ton decays are mediated by k ++ at tree level, and 
set very stringent constraints on the correspond- 
ing Yukawa coupling g and the mass of k ++ . The 
bounds at 90 % C.L. read 
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Note that NSIs are induced by exchanging the 
singly charged Higgs h + , and one can in princi- 
ple tune the mass of k ++ or the scale of g in order 
to suppress its contributions to the LFV decays. 
However, since all the parameters /, g, nih, and mk 
enter in the expression for the neutrino mass ma- 
trix, c.f. Eq. (H]), the constraints from Eqs. (jTTJ) have 
to be included also in an analysis of NSIs, in order 
to obtain the correct parameter space available for 
the model. 

/i + e~ — > /i~e + . The muonium to antimuonium 
conversion through the exchange of fc ++ are well 
bounded experimentally. For the similar reason 
mentioned above, these constraints are mainly on 
mk and g, and the current bound at 90 % C.L. is 
\g ee g* J < 0.2 (m fc /TeV) 2 . 
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• Universality in t a — > l^Vi* decays. The Fermi cou- 
pling constant measured in muon and tau decays 
obtains corrections from the exchange of h + , which 
sets strong constraints on the Yukawa coupling /: 

|/ eAl | 2 < 0.015 (m h /TeV) 2 , 
\\Ut\ 2 - \fer\ 2 \ < 0.05 (m h /TeV) 2 , 

\\fer\ 2 ~ l/e^H < 0.06 (m h /TcV) 2 , 

\\Ur\ 2 - |/ eM | 2 | < 0.06 (m h /TeV) 2 . (18) 

• Rare lepton decays £~ —> Both h + and k ++ 
contribute to LFV photon interactions at one-loop 
level, and the most stringent bound comes from 
/i — > e'y. Neglecting the contributions from the dou- 
bly charged Higgs, we obtain experimental bounds 
at 90 % C.L. 

\r eT Ur{ 2 < 3.4 x 10- 5 (mft/TeV) 4 , 

\f;j, lT \ 2 < 1.7 (m,/TeV) 4 , 

|/ e y eT | 2 < 0.7 (m./TcV) 4 . (19) 

Note that in our numerical analysis, contributions 
from both singly and doubly charged Higgs are in- 
cluded, see, e.g., Ref. [lol ]. 

Besides the bounds above, there also exist other con- 
straints, like the ii—e conversion in nuclei and the anoma- 
lous magnetic moments of the muon, which are relatively 
loose, and hence will not be considered in this work. The 
perturbativity of the model imposes limits on the Yukawa 
couplings / and g, in particular for very massive charged 
scalars. Similarly, the stability of the vacuum requires 
ii <C 47rmin(m fc ,m fl ) Q. 

IV. NUMERICAL RESULTS FOR NSI IN THE 
ZEE-BABU MODEL 

We have performed a full scan of the parameter space 
of the model in order to obtain predictions for NSIs. Fol- 
lowing Refs. [1, [l(| , we take as independent parameters 
the lepton mixing angles, the Dirac and Majorana phases 
5 and a, the Yukawa couplings g ee , g efl , g eT and one of the 
three fap's, and the scalar masses rrih, mk, as well as the 
trilincar coupling ii. (Neutrino mass-squared differences 
are fixed to their best-fit values, since their uncertainties 
are comparably small.) The remaining Yukawa couplings 
f a p and g a /3 are then fixed by Eqs. (|2|), (fT3")) . and (fl~4"l) 
for NH or Eqs. Q, |T5]), and for IH. For each set of 
these parameters, we compare the model predictions to 
the experimental data with a \ 2 function 

x^E^r^' (20) 

where p° represents the data of the ith experimental ob- 
servable, <Tj the corresponding la absolute error, and 



Pi the prediction of the model. The experimental ob- 
servables are the neutrino mixing angles (taken from 
Ref. [13]), and the constraints from LFV and universality 
tests given in Eqs. (fT7 ) -([T9" ]) . 

For the dimensionful parameters m^m^, and /i, we 
adopt first two representative choices, namely mh — 
mk = M = lOTcV or rrih = mk = jx = ITeV. 
The latter case might be just in reach for the LHC. 
For a luminosity of 300 fb _1 and under optimistic as- 
sumptions, this may lead to order 10 four-lepton events 
from the pair-production of the doubly charged scalars 
pp -»■ k++k-- -»■ £+£+£'£- [n|. We do not consider 
much lower scalar masses, since already at 1 TcV many 
of the experimental bounds are saturated. At 10 TeV, 
the constraints are less tight, which leaves more freedom 
in choosing the parameters of the model, with the obvi- 
ous disadvantage of not being testable at LHC. Choosing 
fi of the same order as the scalar masses is conservative 
in the following sense. For fixed scalar masses, neutrino 
masses are proportional to fi. Hence, decreasing ii would 
require to increase the Yukawa couplings, see Eq. (fj|. 
This would make the constraints from Sec. IIIII more se- 
vere and the parameter space would be more constrained. 
Therefore, we decided to take rrih = rrik = // in order to 
keep fj, relatively large, but still ensure the stability of 
the vacuum Q. Towards the end of this section we also 
investigate the case rrih ^ ?^fc- 

In Fig. [21 we present the allowed regions of NSI param- 
eters at lcr, 2a, and 3a C.L., defined as contours in A% 2 
for two degrees of freedom with respect to the \ 2 mini- 
mum. The left (right) panels correspond to scalar masses 
of 10 TeV (1 TeV). The upper plots in Fig. H show the 
NSI parameters |e^J and |e^-| for neutrino propagation. 
We find that values of order e ~ 10~ 3 can be obtained 
only in the IH case, whereas in the NH case they are 
typically one order of magnitude smaller. These flavor 
diagonal propagation NSI parameters induce effectively 
a non-standard matter effect and it has been shown that 
present long-basclinc experiments are not very sensitive 
to these NSI effects [2l|. Even a two-baseline neutrino 
factory would only be sensitive to such NSI parameters at 
the level of e > 10^ 2 [H]. The lower plots in Fig. [5] indi- 
cate that the NSI parameter \e s eT \ relevant at the source of 
a neutrino factory may reach values up to 10~ 3 for both 
hierarchies. In the IH case, the parameters |e™ T | = |e^ r | 
may also be as large as few xlO -3 . For a scalar mass 
scale of 1 TcV, a non-trivial lower bound on the NSI pa- 
rameters of order 10~ 4 is found in the right-hand column 
of Fig. Indeed, for scalar masses in the TeV range, 
the model is rather constrained and the requirement of 
a correct neutrino mass matrix pushes the Yukawa cou- 
plings close to the bounds from Sec. IIIII 0, [13, which in 
turn implies "large" NSIs. The extended region seen in 
the figure comes mainly from the freedom to adjust the 
Dirac CP-violating phase S. 

Source NSIs related to the muon decay in a neutrino 
factory are probed efficiently with a near detector, see 
e.g., Ref. [23], since they lead to the appearance of 
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FIG. 2: The allowed region of NSI parameters at la, 2a, and 3a C.L. in the Zee-Babu model. We take mj, = rrik = /i = 10 TeV 
for the figures in the left-hand side column and — rrik = fi = 1 TeV for the figures in the right-hand side column. 



"wrong" neutrino flavors even at "zero distance" [24| . 
For the cases of interest in the Zee-Babu model, \e* t \ 
and |e^ T |, obviously a tau detector at the near site would 
be useful 0, ES Si SE S3, S3 • The authors of Ref . [26| 
consider as an example a 2 kt OPERA-like near detector 
and find sensitivities for |e| t |, |e„ r | > 7 x 10~ 4 . Note 
that in order to disentangle the effect of e s eT and T , the 
ability to identify the charge of the tau lepton would be 
required. 

The sensitivity to \e™\ for neutrino pro pag ation has 
been discussed for atmospheric neutrinos |29l [30j , the 
OPERA long-baseline experiment [3l| , superbeam exper- 
iments [IH, and a neutrino factory, e.g. in Refs. [22I. I27I 
[H, [13] ■ Typically, the reach is at a few xlO -2 or worse, 
which is not sufficient to probe the parameter range pre- 
dicted by the Zee-Babu model. However, there are two 



reasons why in our case we may expect better prospects 
to observe NSI effects in this channel. First, in many of 
the above mentioned studies the complex phase of e™ T 
has been marginalized over, whereas in the Zee-Babu 
model e™ r is predicted to be real, see the discussion af- 
ter Eqs. (fT5)) and (fTl)]) . Second, Eq. (JHJ) relates source 
and propagation NSIs in this channel. The relevance of 
the phases can be understood from Eq. (35) of Ref. , 
which shows that the relevant leading terms in the sur- 
vival probability P pp are proportional to |e*| sin(0 s ) and 
\e™ T \ cos(4> m ), where (j> s ' m = arg(e*'™). Hence, these 
terms can be set to zero if 4> s and <\> m can be chosen 
independently, but if they are coupled by <\f = ir — 4> m 
following from Eq. at least one of them will always 
be non-zero. Under these special conditions, we estimate 
from the results of Refs. [2(| S3 sensitivities of a neu- 
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FIG. 3: Allowed regions at la, 2a, and 3a C.L. in the plane of sin 2 #13 and |e;! T | (upper panels) or |e^ r | = |eJT T | (lower panels). 
We take mj = = fi = 10 TeV for the figures in the left-hand side column and mj, = = fi = 1 TeV for the figures in the 
right-hand side column. 



trino factory for e™ T in the range of 10~ 3 , even without 
a near detector. In the presence of a near tau-detector, 
a sensitivity for |s™ | = |e*| > 6 x 10 -4 is reached (26j . 

In Fig. [3l we show the correlations between NSI pa- 
rameters and the mixing angle 6* 13 . In the IH case and 
scalar masses at the TeV scale, one obtains a quite strong 
prediction for the mixing angle #13. From Eqs. (|15p and 
(fT6)) follows that \fn T \ is suppressed by S13, whereas a 
correct neutrino mass matrix requires / MT to be of the 
same order as f eft and f eT . Figure [3] indeed shows that 
for the IH case and scalar masses in the TeV range, val- 
ues of #13 close to its present bound are predicted [l(|, 
with a lower bound of sin 2 (9 13 > 10~ 2 . Such a sizable 
lower bound is of particular interest, since it would guar- 
antee a discovery at the forthcoming reactor [35l , l3a ] or 



long-baseline [3 71 ] experiments in the near future [381 ] . 

In the NH case, no lower bound on #13 is obtained. In 
this case, the presence of NSIs may have an impact on 
the search for #13 at a future neutrino factory, especially 
if 6*i3 is relatively small. In particular, e s er may lead to 
v T at the source, which will oscillate to and lead to 
so-called "wrong-sign" muons in the far detector, which 
might be confused with the effect of a tiny #13 [39[ . 

Finally, let us relax the assumption m% = and 
investigate the dependence on the masses of the scalars. 
In Fig. 01 we show the size of NSI parameters by fixing 
one of the two scalars at 1 TeV and varying the mass 
of the other one. In order to guarantee the stability of 
the vacuum Q , while keeping neutrino masses as large as 
possible, we use for the scale of lepton number violation 
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FIG. 4: Allowed regions at la, 2a, and 3a C.L. in the plane of the scalar masses and |e^ r | = \s™ T \- In the left panel, we vary 
the mass of the singly charged scalar and fix the doubly charged one at = 1 TeV, while in the right panel, we vary mt 
and fix Trih — 1 TeV. In both cases, we use /i = min(mh, m*,). 



\i = min(m/i, mfe). The qualitative behaviour of these 
results can be understood from the expression for the 
neutrino mass matrix in Eq. (j4]) and the fact that the 
constraints on the Yukawa couplings / and g from Sec. IIIII 
scale with irih and rrik, respectively. 

First, the lower bound on the scalar masses follows 
from the fact that decreasing either or m& decreases 
the neutrino mass because of fi = mii^m/j, m^). At the 
same time, the bounds on the Yukawas become more se- 
vere and it is impossible to obtain sufficiently large neu- 
trino masses. Second, if we increase rrik, while keeping 
mh at 1 TeV (right panel of Fig. [3J , neutrino masses de- 
crease with m^ 2 because of M = max(m/,,mjfe). Since 
the constraints on g increase only with m^, it is not 
possible to compensate for the m^ 2 suppression by in- 
creasing g. Hence, if the singly charged scalar is at the 
TeV scale also the doubly charged one has to be in that 
range. However, the opposite statement is not true and 
relatively large rrih is possible for = 1 TeV. If m/j is 
increased, again neutrino masses decrease as m^ 2 , but 
in this case there are two factors of / entering Eq. (01, 
and since the constraints on them increase as mh, it is 
possible to keep neutrino masses constant by increasing 
/. However, note that this does not lead to larger NSIs, 
since e oc / 2 / m /l> see Eq. which is constant, in agree- 
ment with the left panel of Fig. [4j 

We conclude that from the point of view of NSIs, 
the crucial mass parameter is the one of the doubly 
charged scalar. Note also that this one has the most 
striking signature at colliders, namely the decay into two 
like-sign leptons. Apparently, vastly separated masses 
for the singly and doubly charged scalars is either not 
phcnomcnologically viable or does not affect the predic- 
tion for NSIs. Therefore, the NSI results obtained for 



mh = rnk are generic. 



V. DISCUSSION AND SUMMARY 

We have studied NSIs in the Zee-Babu two-loop neu- 
trino mass model, which are mediated by the singly 
charged scalar of the model. We have shown that non- 
standard neutrino matter interactions relevant for the 
propagation of neutrinos may be induced. The relevant 
parameters e™ , e™ , and e™ T may reach values of order 
10~ 3 . While flavor diagonal NSIs of this size are too 
small to be observable, the off-diagonal term e™ T may be 
within the reach of a future neutrino factory. In addition 
to these matter NSIs, NSIs affecting the muon decay at 
the source of a neutrino factory may be induced in the 
Zee-Babu model, in both the v e — > v T (e| T ) and — > v T 
(e^ T ) channels. 1 The possible size of NSI parameters de- 
pends on the scale of new physics (i.e., the masses of 
the singly and doubly charged scalars, and m,k, re- 
spectively, and the scale of the lepton-number violating 
parameter /i) and on the type of the neutrino mass hier- 
archy, NH or IH. 

The most constrained situation is obtained for IH and 
scalar masses at 1 TeV. In this case, the NSI parameters 
e s eT and are predicted to be in the range 10~ 4 — 10~ 3 , 
probably in reach of a near tau-detector at a future neu- 
trino factory [26|. Thanks to the fact that e* r is real 



Note that in superbeam or reactor experiments the neutrino pro- 
duction proceeds via hadronic interactions, which are not af- 
fected in this model. 
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in this model combined with the relation e* = — e™ T *, 
even a standard neutrino factory without a near tau- 
dctcctor may be sensitive to the NSI values predicted 
in this case. Furthermore, this configuration predicts a 
value ofj9i3 in reach of the upcoming oscillation exper- 
iments [38[, as well as signals in LFV processes close to 
the present bounds, with good prospects for a signal in 

If kinematically accessible, the singly and doubly 
charged scalars of the model could be directly produced 
through the s-channel processes at the Tevatron and the 
LHC. There is no severe suppression of the cross section 
for the production of doubly charged scalar k ++ , and if 
2m h + > m k ++ , it will predominantly decay into like-sign 
charged-lepton pairs with a very striking experimental 
signature. However, note that doubly charged scalars oc- 
cur in a variety of models (e.g., the triplet scalar model 
for neutrino mass), and therefore, complementary signa- 
tures are required to identify the model. If the singly 
charged and doubly charged scalars are found at LHC, 
then — besides signals in LFV searches — the Zee-Babu 
model predicts a rather large value of #13 and signals 
for NSI at a neutrino factory at a level of 10~ 3 if the 
mass hierarchy is inverted. In case of NH, #13 as well as 
NSIs may still be in reach of future experiments, but no 
signal is guaranteed in cither case. A similar situation 
emerges if the scale of new physics is beyond the reach of 
the LHC, for example at 10 TeV. In this case, observable 



signals may still arise for NSIs at a neutrino factory (as 
well as for #13), but no relevant lower bound is obtained. 

In conclusion, we have shown that, in the case of the 
Zee-Babu model for radiative neutrino masses, the in- 
terplay of the phenomenology at colliders, searches for 
LFV, and NSI effects at a neutrino factory could play a 
complementary role towards the goal of identifying the 
true mechanism of neutrino mass generation. 
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